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We study a massive cosmic strings with BII symmetries cosmological models in two contexts. The 
first of them is the standard one with a barotropic equation of state. In the second one we explore 
the possibility of taking into account variable "constants" (G and A). Both models are studied under 
the self-similar hypothesis. We put special emphasis in calculating the numerical values for the 
equations of state. We find that for u £ (0, 1], G, is a growing time function while A, behaves as 
positive decreasing time function. If ui — 0, both "constants", G and A, behave as true constants. 



o 
o 

(N 



(N 



cr 

i , 



> 
in 

o 

On 
O 



X 



I. INTRODUCTION 

The exponential expansion of the Universe (inflation- 
ary era) causes the Universe to heat up to a very high 
temperature so the subsequent evolution of the Universe 
is exactly as in hot BB model. Hence, the phase transi- 
tion (as the temperature falls below some critical temper- 
ature) in the early universe causes topologically stable de- 
fects: vacuum domain walls, strings and monopoles (see 
[l| and H). But domain walls and monopoles are disas- 
trous for the cosmological models. Strings, on the other 
hand, causes no harm, but can lead to very interesting 
astrophysical consequences (see Q). Also the existence 
of a large scale network of strings the early universe does 
not contradict the present-day observations. The vac- 
uum strings may generate density fluctuations sufficient 
to explain the galaxy formation (see Q). 

The relativistic treatment of strings was initiated by 
Letelier (see @-@) and Stachel (see Here we have 
considered gravitational effects, arised from strings by 
coupling of stress energy of strings to the gravitational 
field. Letelier (see @) defined the massive strings as the 
geometric strings (massless) with particles attached along 
its expansions. 

The strings that form the cloud are massive strings 
instead of geometrical strings. Each massive string is 
formed by a geometrical string with particles attached 
along its extension. Hence, the string that form the cloud 
are the generalization of Takabayasi's relativistic model 
of strings (called p— strings). This is simplest model 
wherein we have particles and strings together. In princi- 
ple we can eliminate the strings and end up with a cloud 
of particles. This a desirable property of a model of a 
string cloud to be used in cosmology since strings are not 
observed at the present time of evolution of the universe 
(see 

In modern cosmological theories, the cosmological con- 
stant remains a focal point of interest (see [HI]- 03 f° r 
reviews of the problem). A wide range of observations 
now compellingly suggest that the universe possesses a 
non-zero cosmological constant. Some of the recent dis- 
cussions on the cosmological constant "problem" and 
on cosmology with a time- varying cosmological constant 
point out that in the absence of any interaction with 
matter or radiation, the cosmological constant remains 



a "constant". However, in the presence of interactions 
with matter or radiation, a solution of Einstein equa- 
tions and the assumed equation of covariant conserva- 
tion of stress-energy with a time- varying A can be found. 
This entails that energy has to be conserved by a de- 
crease in the energy density of the vacuum component 
followed by a corresponding increase in the energy den- 
sity of matter or radiation Recent observations strongly 
favour a significant and a positive value of A with mag- 
nitude A(Gh/c 3 ) w 1CU 123 . These observations suggest 
on accelerating expansion of the universe, q < 0. 

Our current understanding of the physical universe is 
anchored on the analysis of expanding, isotropic and ho- 
mogeneous models with a cosmological constant, and lin- 
ear perturbations thereof. This model successfully ac- 
counts for the late time universe, as is evidenced by the 
observation of large scale cosmic microwave background 
observations. Parameter determination from the anal- 
ysis of CMB fluctuations appears to confirm this pic- 
ture. However, further analyses seem to suggests some 
inconsistency. In particular, its appears that the uni- 
verse could have a preferred direction. Followup anal- 
yses of various sets of WMAP data sets, with different 
techniques seem to lead to the same conclusion. It is 
still unclear whether or not the directional preference 
is intrinsic to the underlying model, and what implica- 
tions this has on our understanding of cosmology. For 
this reason Bianchi models are important in the study of 
anisotropics. 

The study of SS models is quite important since a 
large class of orthogonal spatially homogeneous models 
are asymptotically self-similar at the initial singularity 
and are approximated by exact perfect fluid or vacuum 
self-similar power law models. Exact self-similar power- 
law models can also approximate general Bianchi models 
at intermediate stages of their evolution. This last point 
is of particular importance in relating Bianchi models to 
the real Universe. At the same time, self-similar solu- 
tions can describe the behaviour of Bianchi models at 
late times i.e. as t — > oo (see [l9j]).A particular interest is 
in determining the exact value of the equations of state 
(see Hil). The geometry and physics at different points 
on an integral curve of a homothetic vector field (HVF) 
differ only by a change in the overall length scale and 
in particular any dimensionless scalar will be constant 
along the integral curves. In this sense the existence of a 
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HVF is a weaker condition than the existence of a KVF 
since the geometry and physics are completely unchanged 
along the integral curves of a Killing vector field (KVF) . 
However, the existence of a non-trivial HVF leads to re- 
strictions on the equations of state. 

Therefore the paper is organized as follows. In sec- 
tion II we outline the model as well as the self-similar 
solution. In section III we shall study a massive cosmic 
string model with barotropic equation of state under the 
self-similar hypothesis. We put special emphasis in cal- 
culating the values of the equation of state which made 
that the solution follows a power law. In section IV we 
shall study a massive cosmic string model with barotropic 
equation of state under the self-similar hypothesis but al- 
lowing that "constants" G and A may vary. Section V is 
devoted to study the curvature behaviour of the model 
taking into account the obtained solutions in the above 
sections. We end summarizing the main results. 



II. THE MODEL. 

In synchronous co-ordinates the metric is: 

ds 2 = -c 2 dt 2 + a 2 (t)dx 2 + (b 2 (t) + K 2 z 2 a 2 (t)) dy 2 + 
+ 2Ka 2 (t)zdxdy + d 2 (t)dz 2 , (1) 

where the metric functions a(t), b(t), d(t) are functions of 
the time co-ordinate only and K G R. The introduction 
of this constant is essential since if we set K = 1, as it 
is the usual way, then there is not SS solution for the 
outlined field equations. In this paper we are interested 
only in Bianchi II space-times, hence all metric functions 
are assumed to be different and the dimension of the 
group of isometries acting on the spacelike hypersurfaces 
is three. 

Once we have defined the metric and we know which 
are its killing vectors, then we calculate the four velocity. 
It must verify, C^Ui — 0, so we may define the four 
velocity as follows: 



"' = [ ^0,0,0) 



in such a way that it is verified, g(u % , u 1 ) = —1. 
From the definition of the 4- velocity we find that: 

. i 1 fa' V d'\ 1 
9 = u 1 ., = - - + - + - = -H, 
c V a b d I c 



(2) 



i 
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[H 2 + H 2 + H 2 — H 1 H 2 - H 1 H 3 - H 2 H 3 ) 



From equation CvQij = ^9ij, where V € X(M), we find 
the following homothetic vector field for the BII metric 



b' 



V = td t + 1 - t-\ xd x + 1 - t-\ yd y + 1 - t-\ zd z , 



d' 



with the following constrains for the scale factors: 



a{t) = a t° 



b(t) = b t° 



d{t) = d t° 



(3) 



(4) 



with <2i, d2, 03 € R, in such a way, that the constants at 
must verify the following restriction (see [HI): 

a-2 + a>3 — ai = l. (5) 
Taking into account the field equations (FE) 
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Rij - 2 R 9ij - —^r Ti 3 _ V'/..- 



(6) 



where A is the cosmological constant, and the energy- 
momentum tensor, X^-, for a cloud of massive strings is 
given by 



Tj = (p + p) u J u t + pgf - Xx J Xi, 



(7) 



where p(t) is the rest energy density, i.e. is the rest en- 
ergy density of the cloud of strings with particles attached 
to them (p— strings). X(t) is the string tension density, 
which may be positive or negative, u 1 is the four-velocity 
for the cloud particles. x % is the four-vector which rep- 
resents the strings direction which is the direction of 
anisotropy and p — p p + A, where p p denotes the par- 
ticle energy density (is the cloud rest energy density), 
i.e. the string tension density is connected to the rest 
energy p for a cloud of strings (p— strings) with particle 
attached to them by this relation (see [9| and flO]). 

Since there is no direct evidence of strings in the 
present-day universe, we are in general, interested in con- 
structing models of a Universe that evolves purely from 
the era dominated by either geometric string or massive 
strings and ends up in a particle dominated era with or 
without remnants of strings. 

Moreover the direction of strings satisfy the standard 
relations: 



l. 



u l Xi = 0, 
x l = (0,a-\0,0) 

The Raychaudhuri equation reads 
f 



(8) 
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2a z + RijU l u\ 



with 



(9) 



(10) 



It is customary to assume a relation between p and A 
in accordance with the state equation for strings. The 



3 



simplest one is a proportionality relation p = aX, where 
the most usual choices of the constant a are the following 
ones (see [13]). Geometric strings. Nambu. a = 1, p = A, 
Pp = 0. p— strings or Takabayasi strings 



a = l + W, 



(11) 



with W > 0, iff p p — WX. Reddy strings. a — — 1. A more 
general "barotropic" equation of state is p = p (A) , p p = 
p—X. In the case of Takabayasi strings if W is very small, 
then only geometric strings appear. On the other hand, if 
W is infinitely large, then particles dominate over strings. 



III. BAROTROPIC EQUATION OF STATE 

We shall consider the Takabayasi's equation of state, 
i.e. p = aX, with a = 1 + W, so the resulting FE are: 
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(up - A) 
(13) 

up, (14) 
up, (15) 
(16) 



noting that b = d. So once again the model collapses to 
a LRS BII model (see [H]) where 

a(t)=a t ai , b(t) = b t a ' 2 = d(t), 

and the constrain 02+03 — 01 = 1, reads now 2a2 — 01 = 1. 
Therefore the FE yield: 
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(up-X), (18) 
u>p, 
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(19) 
(20) 



In this way, following the standard procedure (see [231 ]) 
we find that 

9 = Pot- f3 , (3 = (4o 2 - 1) (a, + 1) - ^jp-, 



and hence we need only to solve the system of algebraic 
equations 



3K 2 



3a2 — 202 — = —A [ LO — 



1 



7a2 - 8o 2 + 2 + ^- = -uA, 



K 2 
4 



(4a a -l)( W + l)-^^ = L>. 



1 + W 



1 + w J ' 

(21) 

(22) 
(23) 



with 



A E 2 o ^ ^ -° 

A = 5o, — 2a 9 . p n = — ^, A= . 

2 2 4 ' P0 8ttG' 1 + W 



By solving the system (|2HI23[) we obtain the following 
results: 



o 2 



2 + co(l + W) + 3W 
2{l + 2W + 2uj{1 + W))' 



(24) 



and 



K 



(W - 2) 2 + 2uj (W 2 + 3W + 2) - 3uj 2 (1 + W) 2 - 



[1 + 2W + 2u{l + W)) 



(25) 

We may list the physical solutions for different values 

Of LO 



1. If we set u) = 0, then we get: 



yiy 2 -4(iy + i) 
x = 1T2M 7 ' a2 



2 + 31^ 
2(1 + 2VF)' 



and the following restriction: W > 2 + 2%/2. There- 
fore K e (0,0.5), while o 2 € (0.75,0.77) and 
oi € (0.5,0.54). 



2. If w = 1/3, 



K = 



V-27 - 24IT + 12W 2 



K G (0,0.43301), 



with W > 1 + ^ = 2.8028, therefore the scale 
factor 02 behaves as: 

7+ 10VK 



"2 



o 2 G (0.625,0.638), 



2 (5 + 8W) 
therefore 01 G (0.25, 0.277) . 

In the following table we summarize the behavior of 
the main quantities for this model 
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W 


K(W) 


a 2 (W) 


ai(W) 


OBS 


1 


< 








unphysical 


1/3 


> 2.8028 


(0,0.43301) 


(0.625,0.638) 


(0.25,0.277) 


physical 





> 4.8284 


(0,0.5) 


(0.75,0.77) 


(0.5,0.54) 


physical 


-1/3 


< 








unphysical 


-2/3 


< 








unphysical 



and the energy density behaves as: 

P = Pot~^ /3=(4 02 -l)( £l ; + l)-^=il«2. 

Therefore we are only able to obtain solutions for lu = 
1/3 and u> = 0. For the rest of equations of state the 
obtained solution is unphysical since in all these cases we 
have obtained W < 0, in order to get K € K. But we 
have ruled out the equation of state W < 0. 

With the obtained results we can see that 



H = 



4a 2 - 1 



t 



2- 4a 2 
q = < 0, 



4a 2 



1 



note that 2 — 4a 2 < 0, Vw e [0,1/3] , while the shear 
behaves as: 



(ga - If 
3c 2 t 2 



0. 



As it is observed, these quantities fit perfectly with 
the current observations obtained by High -Z Supernova 
Team and Supernova Cosmological Project (see for ex- 
ample [II, [iff, [11 and (H). 

IV. VARIABLE CONSTANTS. 

In this model the resulting FE are: 
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(26) 

Ac 2 , 
(27) 
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noting that b = d. So once again the model collapses to 
a LRS BII model (see [Hj]), where 

a(i) = a t ai , b(t) = M" 2 = 

and the constrain 02+03 — 01 = 1, reads now 2a2 — 01 = 1, 
as above. Therefore the FE yield: 
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(wp - A) + Ac 2 , 
(33) 

(34) 
(35) 
(36) 



From eqs. (f3"5j) and ([57]) we get 

P = Pvt-\ 7=(4aa-l)(w+l)- 



(2o 2 - 1) 
1 + W 



(37) 



(38) 



where 7 = (w + 1) a — y^fy and a = (ai + 202), 2a2 
01 = 1, and we shall assume that po > 0. 
From eq. (f3"2"]) we obtain: 



A =- 

c z 



9 87rG 
At~ 2 ^Pvt 1 



(39) 



where A = 2aia2 + a 2 — 

Now, taking into account eq. ([56")) and eq. (|3"9")) . algebra 
brings us to obtain 



G = Gnf ~ 2 , 



Gn — 



c 2 A 



47T/907 ' 

and therefore the cosmological "constant" behaves as: 



(40) 



(41) 
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With all these results, we find that the system to solve 
is the following one: 



3a 2 — 2a 2 



3K 2 2A 
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1 



4 — r + 1) "2 1 + W 



(42) 

7a2-8a 2 + 2+^ = -^(( W + l)-|), (43) 

and whose solutions may be listed as follows: 
1. If uj = 0. 

oi = 2a 2 - 1, (44) 

and 



_ ^(1^-4(^ + 1)) /_ 1 \ 
W(1 + 2W) fc ^ 2'7 ' 

therefore this solution has only sense if W > 2 + 
2V2 « 4.8284. In this way we find that 

and the behaviour for the rest quantities is: 

G«i 7 ~ 2 = G , G = = const. > 0, 

47rp 7 

A«i" 2 , A = 4fl-- N ) =0, 



with A = a 2 (5a 2 - 2) - ^ > 0. Note that this 
solution coincides with the above one where the 
constants behaves as true constants and A vanishes. 



2. If uj = 1 



K = - 



V-27W 2 - 22W - 7 - 12VK 3 + AW 4 



{-1 + W) (3 + AW) 
so K £ (-5,0) with > 4.7, and 
-1 + 3I^ + 61^ 2 



02 



2(-W-3 + 4W 2 



e (0.75,0.90), 



hence ai £ (0.5, 0.80) , with A = a 2 (5a 2 - 2) - 
^ « 1.25 > 0, and therefore 

p~t\ 7 = 2(4a 2 - 1) - {2 ^~^ €(4,5.06), 



G » ^" 2 , Go 



c 2 A 1 



47T/90 7 



>0, 



3. If uj = 1/3, we get the following solution: 

V-519VK 2 - 78W - 15 - 456W 3 + 144VK 4 



if = -- 



(-1 + W) (5 + 8VK) 



so K £ (-5,0) , with > 4.1, and 
-l + 23iy + 36iy 2 



a 2 



G (0.75,0.81). 



2 (7W - 5 + 24W 2 ) 
hence a\ £ (0.5,0.63), therefore 

p«t\ 7 = \ (4a 2 - 1) - {2 * 2 + -^ ] £ (2.66,2.86) : 



G w F- 2 , Go = 



c 2 A 1 



47rp 7 



>0, 



A = - 1 - - > 



\ 7, 

4. If w = -1/3, then 

V-243I^ 2 -6W + 9 - 228W 3 + 36W 4 



K 



(1 + 3W)(1 + 4W) 
so K £ (-i, 0) with W > 6.6, and 
-1 + 11VK+18T^ 2 



a 2 



e (0.75,0.7507), 



2 (7W + 1 + 12W 2 ) 
hence ai € (0.5,0.501), therefore 

p»f, 7= g(4a 2 -l)- ^^ g (1.26,1.33), 



G » F- 2 , Go 



c 2 A 1 



47rpo 7 



>0, 



A«^ 2 , A = ^(l--]<0. 



A / 2 

c 2 \ 7 



5. If u = -2/3, 

V-1321^ 2 - 48W + 48 - 132VK 3 + 9W 4 



K = — 



{2 + 3W) (-1 + 2W) 



in such a way that K £ (-|,0) with > 15.2, 
and 

-A + 2W + 9W 2 
°^ 2(^-2 + 6^) e '°- 75 ' a752 '' 

hence 01 e (0.5,0.504), therefore 

7= I (4a 2 -l)-^-^e (0.63,0.66), 



G w T" 2 , Go 



c 2 A 1 



47rpo 7 



>0, 



A * / 7 A„ = 4 fl- ^ ) <>• 



In the following table we have summarized the behav- 
ior of all the quantities 
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K{W) 


a 2 
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G( 7 -2) 
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OBS 
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> 4.7 


(-0.5,0) 


(0.75,0.90) 


(0.5,0.80) 


(4,5.06) 


(2,3.06)/ 


>o\ 


phys 


1/3 


> 4.1 


(-0.5,0) 


(0.75,0.81) 


(0.5,0.63) 


(2.6,2.8) 


(0.6,0.8)/ 


>o\ 


phys 





> 4.8 


(-0.5,0) 


(0.75,0.77) 


(0.5,0.54) 


2 


const. 





phys 


-1/3 


> 6.6 


(-0.5,0) 


(0.75,0.7507) 


(0.5,0.501) 


(1.26,1.33) 


(-0.74,-0.67) \ 


<o\ 


phys 


-2/3 


> 15.2 


(-0.5,0) 


(0.75,0.752) 


(0.5,0.504) 


(0.63,0.66) 


(-1. 37,-1. 34) \ 


<o\ 


phys 



Therefore, as we can see, and comparing this table 
with the above one obtained in the last section we may 
conclude that to allow that the "constants" may vary en- 
large the set of solutions. In this case we have obtained 
a solution for equation of state ui £ (—1,1] while in the 
standard solution we only have a solution for u> = 0, and 
lu = 1/3. In the same way it is interesting to empha- 
size that for tu > i.e. to £ (0, 1], the "constant", G, is a 
growing time function while the cosmological "constant" , 
A, behaves as positive decreasing time function (see for 
example [l5[ , [Ty] , [13] and [HI). For lv = 0, we have that 
both "constants" , G and A, behave as true constants, in 
fact we have obtained the same solution as in the above 
section. If to < i.e. u> £ (—1,0), the "constant", G, is 
a decreasing time function while the cosmological "con- 
stant" , A, behaves as negative decreasing time function, 
so, from the observational data we may rule out these 
solutions. 

With the obtained results we can see that 



H = 



4a 2 



1 



2-4a 2 



t 



4a 2 



1 



<0, 



note that 2 - 4a 2 < 0, and 4a 2 - 1 > 0, Vcj £ [0, 1] , while 
the shear behaves as: 



(« 2 -ir 

3c 2 t 2 



0. 



As it is observed, these quantities fit perfectly with 
the current observations obtained by High -Z Supernova 
Team and Supernova Cosmological Project (see for ex- 
ample 0,(3, [HI and 0). 



V. CURVATURE BEHAVIOUR. 

With all these results, we find the followin g be haviour 
for the curvature invariants (see for example [24|-p7|). 
Ricci Scalar, I , yields 



2 / K 2 c 2 
/ o = ?72 (llal,-10a 2 + 2- — 

while Krestchmann scalar, I\ := R^kiR 1 ^ 1 , yields: 
1 



(45) 



h = T -^ J (432a 4 - 960fl3 -(- 896a 2 - 384a 2 + 64 

2 



4c 4 < 4 

+K 2 c 2 (UK 2 c 2 - 40a^ + 80a 2 - 48)) 



The full contraction of the Ricci tensor, 7 2 := 



h = 



1 



(528a^ - 1024a^ + 768a 2 2 - 256a 2 + 32 



4c 4 i 4 

K 2 c 2 (3ifV-16a 2 + 8)) , 

this means that the model is singular. 

The non-zero components of the Weyl tensor. The 
following components of the Weyl tensor run to ±oo when 

t -> 0, 



and 



r< ~ r< ~ +4(a 2 -l) /-* ^ i 2 (a 2 -l) 
^txtx ~ <^txty ~ 1 ) ^tztz ~ I i 



these others run to zero when t — * 0, 

(~* ~ r 1 ~ r' ~ j.2(3a 2 — 2) 

^xzyz ~ ^xzxz ~ ^xyxy ~ 1 ; 



(-1 ^ f 2(2a 2 -l) , 2.2(3a 2 -2) 

^yzyz ~ L t ^ <- 



txyz ~ ^tyxz ~ ^tyyz ~ ^tzxy ~ ^ ) 



J txyz ~ ^ty 

noting that when a 2 — > 3/4, then they run to a const 
~ I 

3 



The Weyl scalar, 7 3 = C abcd C abcd = h - 2/ 2 + \l 2 , 



(this definition is only valid when n = 4) 

^3 = ([(as - 1) (-2 (oa - 1) - 3Kc) + K 2 c 2 } 

x [(o a - 1) (-2 (oa - 1) + 3Kc) + K 2 c 2 ] ) . 

The electric part scalar I4 = EijE^ , (see W.C. Lim et 
al [H) 



I A = 



1 



-2(a 2 -l) +#V , 



6c 4 i 4 

while the magnetic part scalar J5 = HijH^ , yields 

3 ' -\2 .X2 



2c 2 t 4 



(aa-r/i^. 



The Weyl parameter (see W.C. Lim et al (28[) which 
is a dimensionless measure of the Weyl curvature tensor 



W 2 1 1 

W 2 = — = — (/•;,,// • II, ill' 1 ) = ^TTTrC a6c<1 C Q , 



ff 4 6H 4 



2AH 4 



bcd-i 



7 



where 



W 1 



l F ((a 2 -l) 2 + ^ 2 c 2 ) 2 ((2 a2 -2) 2 + ifV) 2 , 



36c 4 i 4 



note that the value of W 2 is really small. We may calcu- 
late the quantity 



W 1 



((aa - l) 2 + K 2 c 2 ) 2 ((2oa - 2) 2 + if 2 c 2 ) ' 
36c 4 (4a 2 - l) 4 



const, 



but as it is observed W — * 0. W can be regarded as 
describing the intrinsic anisotropy in the gravitational 
field. Hence we may conclude that the model is close to 
be isotropic, (sec [20]) since we have that it is verified 
the Weyl isotropization criterion i.e. W — > 0. 
The gravitational entropy (see [25|-[26]) 



I n 



p A = — = 
/ 2 



/i - 2J 2 + i/ 6 



A 1 4 

= — + r — - 2 = const 

1 2 3 I 2 



since we are working with a SS solution (see [3(| and 
[3l| for a discussion), and all the dimensionless quantities 
remain constant along timelike homothetic trajectories as 
the Weyl parameter. 



VI. CONCLUSIONS 

We have studied two massive cosmic string Bianchi 
type II cosmological models under the self-similar hy- 



pothesis. In the first of the studied models (those with 
"constant" constants) we have obtained that the self- 
similar solution is only valid if the equation of state is 
us = and us — 1/3. For the rest of possible values of 
us the solution is unphysical. In the same way we have 
shown that such solutions are only valid if W (the param- 
eter in the equation of state for the strings) is W > 2.8 (if 
us = 1/3), and W > 4.8 (if us = 0). So may conclude that 
the self-similar solution is quite restrictive. Furthermore 
the obtained solution collapses to a LRS BII solution 
since two of the scale factors are equal. 



In the second of the studied models (those with "vary- 
ing" constants) we have shown that considering a time 
varying constants we may enlarge the set of solutions. In 
this case we have obtained a mathematical solution for 
us G (—1,1], and with similar restrictions for W. In the 
same way it is interesting to emphasize that for us > i.e. 
us E (0, 1], the "constant", G, is a growing time function 
while the cosmological "constant" , A, behaves as positive 
decreasing time function. This fact fits perfectly with the 
current observational data. For us = 0, we have that both 
"constants", G and A, behave as true constants, in fact 
we have obtained the same solution as in the first of the 
studied models. If us < i.e. us e (—1,0), the "constant", 
G, is a decreasing time function while the cosmological 
"constant" , A, behaves as negative decreasing time func- 
tion, so, from the observational data we may rule out 
these solutions. 
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